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We discuss the main features of quantum integrable models taking the classes of universality of 
the Ising model and the repulsive Lieb-Liniger model as paradigmatic examples. We address the 
breaking of integrability by means of two approaches, the Form Factor Perturbation Theory and 
semiclassical methods. Each of them has its own advantage. Using the first approach, one can relate 
the confinement phenomena of topological excitations to the semi-locality of the operator which 
breaks integrability. Using the second approach, one can control the bound states which arise in 
each phase of the theory and predict that their number cannot be more than two. 

I. INTRODUCTION 

Important developments of experimental tools for studying the behavior of quantum matter in low di- 
mensions have triggered a renewed interest in integrable and non-integrable systems. Here I would like 
to mention, in particular, two recent experiments: 

• the experiment on onc-dimensional quantum Ising model, realized by means of CoNb206 (cobalt 
niobatc) [ll and used to probe the Eg symmetry of both the excitations and the correlation 
functions [3, 3 of the system. Such a system was also used to study the phenomena of confinement 
of topological excitations [l-Q which occur once the degeneracy of the vacua is broken by the non- 
integrable deformation of the magnetic field. The quantum Ising model (and its tricritical version) 
will be in the following our paradigmatic models for discussing some of the theoretical predictions 
coming from quantum integrability and its breaking. 

• the experiments on one-dimensional Lieb-Liniger gas, realized by using an optical trap in cold boson 
gas [8l-ll2l|. In the paper [l^l the authors reported, in particular, the preparation of out-of-equilibrium 
arrays of trapped one-dimensional (ID) Bose gases, each containing from 40 to 250 Rb^^ atoms, 
which did not noticeably equilibrate even after thousands of collisions. The result of this experiment 
has recently triggered a lot of interest on the role played by integrability in the off-equilibrium 
dynamics of a quantum system (for a review on the subject, see [l3| and references therein). Notice 
that, although the Bethc Ansatz solution of the Lieb-Liniger model is now almost fifty year old 

[isj , the computation of its correlation functions at equilibrium has remained for long time an 
open problem and only recently a general approach has been set up for computing such quantities at 
zero and at a finite temperature [17| . The Lieb-Liniger model will be our paradigmatic example to 
illustrate the efficiency of quantum field theory methods in the computation of correlation functions 
even in non-rclativistic integrable models. 

The best way to characterize a quantum integrable system is through the non-diffractive scattering of 
its quasi-particle excitations [l^. This approach, closely related to the continuum limit of lattice models, 
directly leads to the formalism of integrable quantum field theory [2^ . It is worth to point out that the 
adoption of the continuum formalism of field theory (integrable or not-integrable) is not only extremely 
advantageous from a mathematical point of view but, in the scaling region nearby the critical points, 
is also perfectly justified from a physical point of view: in such a region the correlation length is larger 
than any other microscopic scale and therefore a universal behavior is expected to emerge. It is such 
a universal behavior that is caught by the field theory. Moreover, a quantum field theory embodies a 
strong set of constraints coming from the compatibility of quantum mechanics with special relativity: 
this turns into general relations, such as the completeness of the multiparticlc states or the unitarity of 
their scattering processes, which are extremely useful to successfully pin down the computation, say, of 
exact matrix elements of the order parameters and the universal shape of the correlation functions. 

The layout of the paper is as follows. In Section 2 we briefiy recall the main properties of quantum 
integrable systems, using the class of universality of the Ising model and the Lieb-Liniger models as basic 
examples of the formalism. Non-integrable models will be the subject of the remaining sections and we will 
present two approaches which have been developed in recent years to study the phenomena accompanying 
the breaking of integrability, such as the confinement of the topological excitations, the decay of particles 
or the nature of the energy spectrum around different vacua. In Section 3 wc will discuss some features of 
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non-integrable models. In Section 4 we present the basic formulas of Form Factor Perturbation Theory, 
based on the exact matrix elements of the integrable models; in Section 5 we will present the main results 
coming from the Semiclassical Method, in particular those relative to the spectrum of neutral excitations 
present in each phase of the system. Our conclusions can be found in Section 6. 



II. QUANTUM INTEGRABLE MODELS 



Our prototype models for the discussion on the main features of quantum integrable systems are the 
Ising model and the Lieb-Liniger model (for recent reviews of these models see 21| and [2J] respectively 
and references therein). As it will become clear below, the key quantity for the solution of these models 
is their S-matrix for this quantity permits to identify the spectrum of their excitations, to compute their 
correlation functions and to determine their thermodynamics properties. 



A. Ising model. 



Consider the Hamiltonian of the one-dimensional quantum spin chain of the Ising model 

i i i 

where cr" arc the Pauli matrices. This model is well known (23 | to have a phase transition at ft. = and Ac = 
1, the latter quantity related, in the corresponding 2d classical Ising model, with the critical temperature 
Tc- At criticality, the Ising model is described by the first unitary minimal models of Conformal Field 
Theory (CFT) [23|- CFT permits to organize the operator content of the continuum limit of the model 
and to set the correspondence 

<->e(r) , af^a(r) (2) 

where e(r) and cr(r) are the primary fields of CFT, identified respectively with the energy and magnetiza- 
tion densities. Hence, in the scaling region nearby the critical point, the Hamiltonian ([1} can be written 
as 

H = (E<'^m+E<) +(^-i)E< + ^E'^^' (3) 

\ i i / i i 

~ HcFT + T J dre{r) + h J dra{r) 
where t = (A — 1) is the displacement from the critical value. 



1. Integrable deformations 

The final expression of the Hamiltonian ([3] ) is particularly convenient for discussing the nature of the 
various theories obtained by varying the coupling constants. Using the null-vector structure of CFT, it 
is possible to show that one gets an integrable model only in two cases Q : 

T 7^ and ft = 0, denoted hereafter as thermal deformation] 
T = and ft 7^ 0, denoted hereafter as magnetic deformation. 

When both couplings are different from zero, the null-vector structure of CFT excludes instead the 
existence of any conserved charges and therefore the system is non-integrable [26j . The integrability of 
the two cases above means that in both thermal and magnetic deformation there is an infinite number 
of conserved charge Qs of spin s, although the charges and their spin are different in the two cases: in 
particular, for the (infinite) spectrum of the spins of the conserved charges we have 



5 = 1, (mod 2) for the thermal deformation; 

s = 1, 7, 11, 13, 17, 19, 23, 29, (mod 30) for the magnetic deformation. 



(4) 
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The first values of tlic spins can be identified with the Coxeter exponents of the algebra SU (2) and Es 
respectively: their appearance is related to the possibility of having alternative coset constructions of the 
Ising model at the critical point (for this point see, for instance, [20|). 

The existence of an infinite set of conserved charges ensures the elasticity and the factorization of 
the scattering processes of the excitations. Moreover, these conserved charges strongly constraint the 
spectrum of the excitations, whose number and mass can be obtained by closing the bootstrap of the 
S'-matrix 0. This approach consists of the following steps: 

• parameterizing the dispersion relations of the particle a by Ea = rua cosh 6a , Pa = ni-a sinh 9a 
(where 9 is the rapidity), the two-body elastic scattering amplitude Sab{Sab) relative to the process 
a X b ^ a X b depends, for the relativistic invariance, on the difference of the rapidity 9ai, = da — Of, 
and satisfies the unitary and crossing relations [60j 

Sab{e)Sab{-0) = 1 , Sab{e) = Sab{i^ - 9) . (5) 



• Poles in the amplitude Sab{0) signal the presence of bound states, either in the s or in i channels. If 
9ab = iu^b identified with the s-channcl pole, the mass of the corresponding bound state is given 

by 

— "^a + ™& + '2-mamb cos wjjf, . (6) 



• The bootstrap principle asserts that bound states are on the same footing of the asymptotic states. 
This principle has a series of important consequences. Firstly, if there is the pole of the particle 
c in the scattering amplitude Sab{9), then there will be the pole relative to the particle a in the 
scattering amplitude Sbc{d), as well as there will be the pole relative to the particle b in Sac{9)- 
Secondly, the position of the poles in the three channels satisfies 

<b + <c + = 2^ . (7) 
and, finally, the amplitudes are related by the bootstrap equation (uJJ^^ = tt — 

Sab{0) - Sac{e + lu'i,)Sad{0 - ZugJ . (8) 



So, to summarize the rules of the game, in order to determine the S-matrix (and the spectrum) of 
an integrable theory by means of the bootstrap approach one has to find a set of poles relative to all 
amplitudes Sab which can be interpreted in terms of bound states of the asymptotic particles and which 
are compatible with the bootstrap equation ([5]). The masses of the particles are determined by the relation 
([5]). In practice all this means starting from the amplitude that involves the lighest particle, therefore 
with the simplest pole structure, and then iteratively applying the bootstrap equations ([8]) to get the 
scattering amplitudes involving the bound states of higher mass. Let's sec how this formalism works for 
the two integrable deformations of the Ising model. 

Thermal deformation. The thermal case is known to correspond to a theory of a free Majorana 
fcrmion of mass m (proportional to t), with two-body S'-matrix simply given by S ~ —1. Such a theory 
does not have additional bound states. 

Magnetic deformation. The magnetic case is more subtle. The set of conserved charges implies the 
existence of a lowest mass particle Ai (which is bound state of itself), plus at least two additional particles 
A2 and ^3; moreover, also thanks to the conserved charges, it is also possible to pin down the location 
of their poles in the amplitude Sii{9). With all these data, the amplitude Sii{9) is given by 



12 3 



^11 (^) - 




^ tanh \ {9 + iirx) 
tanh ■^{9 — inx) 

The pole at 9 ^ corresponds to the particle Ai itself, 9 = corresponds to the particle A2 whereas 
the one at 6' = ij^ to the particle A3. If rrii is the mass of the particle Ai, the mass of these new particles 
are TO2 = 2mi cos ^ and — 2mi cos respectively. 
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Applying the boostrap equation relative to the poles, one could get the amplitudes (6*) and 5*13(6'), 
where other poles appear in correspondence to new particles. Applying the bootstrap equations also to 
these poles and carrying on all the necessary steps to close consistently the process (i.e. a consistent 
interpretation of all pole structure of the amplitudes in terms of the identified particles), A. Zamolod- 
chikov showed that the bootstrap procedure closes with 36 amplitudes and 8 particles, whose exact mass 
spectrum is given by [2| 



mi 
m2 


= m 

^ 2mi cos 


TT 

5 " 


= (1.6180339887..) mi 


"^3 


~ 2mi cos 


TT 

30 


= (1.9890437907..) mi 


m4 


= 2m2 cos 


7n 
30 


= (2.4048671724..) mi 


7715 


= 2m2 cos 


2tt 
15 


^ (2.9562952015..) mi 


mg 


= 2m2 cos 


TT 

30 


= (3.2183404585..) mi 


TO7 


~ 4m2 cos 


lcos— = (3.8911568233..) TOi 
5 30 


ms 


— 4m2 cos 


TT StT 

-cos — = (4.7833861168..) TOi 
5 15 



(10) 



Notice that the mass of the particles A4, . . . , are above the mass threshold 2m: in a generic theory 
these particles would be all unstable whereas here their stability is ensured by the integrability of the 
theory. 

2. Form Factors and Correlation Functions 

In addition to the exact expression of the 5- matrix, in quantum integrable theories one can also compute 
the exact matrix elements of the order parameters on the asymptotic state (the so-called Form Factors) 

F^,,a.,....aAOa,,Oa.,...,e,J = (0 1 $ (0) | (0„ J A,, (0, J . . . (0, J) . (11) 

The Form Factors satisfy a set of functional and recursive equations based on unitarity, crossing, fac- 
torization and pole structure of the 5-matrix: different solutions of these equations identify the operator 
content of the theory [2^ . We refer the reader either to the original references [13, for the key equa- 
tions of the Form Factor approach or to the book (20j for a general introduction to this formalism. Below 
we will simply quote the main results relative to the two integrable deformations of the Ising model. No- 
tice that, once we have determined the form factors of a given operator, its correlation functions can be 
written in terms of the spectral representation series using the completeness relation of the multi-particle 
states. For instance, the two-point correlation function of the operator 0{x) (in the euclidean space) can 
be written as [6l[ 

{o{x)om - £ / "^^foy" <o\o{x)\e,,...,e„ >in i„ < 0i,...,0„|o(o)|o > 

= f:J^^^!(^l^"(^---^'OPexp(^-m.^ (12) 

where r is the radial distance r — ^ x\ + x\ (Figura[T]). Similar expressions, although more complicated, 
hold for the n-point correlation functions. It is worth stressing some advantages of this method. 

• The integrals present in the spectral series are all convergent, in sharp contrast with the formalism 
based on the Feynman diagrams, in which one encounters the divergences of the various perturbative 
terms. The reason is that the form factors employ ah initio all the physical parameters of the theory 
and therefore the divergences of the perturbative series (which, on the contrary, refers to the bare 
parameters) are absent. 



FIG. 1: Spectral representation of the two-point correlation functions. 



• If the correlation functions do not have ultraviolet singularities particularly violent (this is the 
case of the correlation functions of the relevant fields), the spectral series has an extremely fast 
convergent behavior for all values of mr. In the infrared region, i.e. for large values of mr, this is 
pretty evident from the nature of the series, because its natural parameter of expansion is e~"^^. 
The fast convergent behavior also in the ultraviolet region mr — >■ has instead a twofold reson: the 
peculiar behavior of the n-particle phase space in two-dimensional theories, further enhanced by 
the form factors behavior, as clarified in (20l [soj. 

Form Factors of the thermal deformation. The field theory associated to the thermal deformation 
of the Ising model is known to be self-dual for the exchange r — >■ — r: the energy operator changes 
sign under this transformation, e{x) — > -~e{x), whereas the order parameter a{x) goes into the disorder 
operator ^(x), i.e. a ^ ^. For r > 0, <y{x) is a Z2 odd operator with non-zero fermion number while 
^jl{x) is a Z2 even operator with zero fermion number: hence, <y{x) has non-zero matrix elements only on 
odd number of fermions, while /i(x) has non-zero matrix elements only on even number of fermions. For 
T < 0, the role of the two operators is swapped and the situation is reversed: this symmetry is due to the 
self-duality of the model. For r > 0, the matrix elements of the various operators (between the vacuum 
state |0) and the asymptotic states of n fermions identified by their rapidities 1 01,^2, • • ■ ,^n)) are given 
by 



(0|e(0,0)| 



im sinh for n = 2 

otherwise 



2n+l 



I a(0, 0) . . . , e2n+i) = (»)'"+' n t^^l^ ■ 



1<J 



(O|//(O,O)|0i,...,02„) = {if'^JliB.nYi- 



l<3 

Inserted into the spectral representation, the energy form factors give the (euclidean) correlation function 
(e(a;)e(0)) = m'^{Ki{mx) — K'l{mx)) (where Ki{y) are the modified Bessel functions). The magnetiza- 
tion form factors give rise instead to the Painleve' equation satisfied by the correlators {a{x)a{0)) and 
{^{x)fj,{0)) [12] ■ Notice that the matrix elements of the magnetization operators have a pole at Pi^Pj = in. 
This observation will be important when we discuss the Form Factor Perturbation Theory relative to the 
non-integrable deformation of the thermal theory. 

Form Factors of the magnetization deformation. As mentioned above, the Ising model in a magnetic 
field has quite a rich iS'-matrix, with amplitudes that have poles up to 12th order. In addition to the 
functional and recursive equations, the form factors of this theory also satisfy other recursive equations 
related to the higher poles of the S'-matrix: the relative formulas and the final exact expressions of the 
Form Factors can be found in the papers [1, 0]. Here we only report the exact result relative to the 
one-particle amplitudes Z„ of the spin-spin correlation function 

{aiq)ai-q)) = , ^ + „ , + • • • + , ^ +■■■ (13) 

They are given by the square of the one-particle Form Factors 

Z, = \{0\a{0)\A,{e))\' . (14) 



FIG. 2: Quasi-particle weights of the spm-spm correlation function normalized to the weight of the particle with 
the lowest mass. 



Their actual value depends on the normalization of the operator <7{x) but universal numbers can be 
extracted in terms of the ratios Zk = Z^/Zi. Their values is in Table 1 and their plot is given in Figure 
2. Figure 2 has to be compared with the experimental data reported in Figure 4E of the paper [l[ and 
one could easily see the two plots agree one to the other. 



Zi 


1 


Z2 


0.279242 


Zi 


0.0847955 


Za 


0.0496251 


^5 


0.00885986 


^6 


0.00458335 


Z-j 


0.000656846 


Zs 


0.0000224409 



TABLE I: Universal ratios made of one-particle amplitudes in spin-spin correlator. 



B. Lieb-Liniger model 

Our second example of integrable model is the Lieb-Liniger (LL) model that describes the low- 
temperature properties of one-dimensional interacting non-relativistic Bose gases. Its accurate experi- 
mental realization has stimulated a lot of interest and has opened new perspectives in the field of 
strongly correlated systems. Thanks to the highly controllable experimental set-ups, it is possible nowa- 
days to thoroughly investigate questions of general nature concerning quantum extended systems, such 
as the dynamics of integrable systems in the presence of small non- integrable perturbations (e.g., three- 
body interactions and/or a weak external trapping potential), the issue of thermalization in quantum 
integrable and non-integrable systems and the behavior of various susceptibilities and response functions. 



Following the recent papers jl7| , we will briefly discuss how one can compute the correlation functions 
of this model (at finite density and at finite temperature) using its relation to the relativistic integrable 
Sh-Gordon model. The LL Hamiltonian for N interacting bosons of mass m in one dimension is 

^ = -^E9^+2AE'5(-^-^.)' (15) 

i—l * i<Cj 

and the corresponding (non-relativistic) field theory is given by the quantum non-linear Schrodinger 
model [ill, that employs the complex field ^ and the Lagrangian 



/: = -2;^W + .^hA^^-^V')-A|^r. (16) 



The effective coupling constant of the LL model is given by the dimcnsionless parameter 7 = 2mX/h^n, 
where A > is the coupling entering the Hamiltonian (fT5|) while n ~ N/ L is the density of the gas {L is 
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the length of the system). It is convenient to express the values of temperature in units of the temperature 
T]j = h^n^ /2mkB of the quantum degeneracy, r = T /T^,. The non-relativistic two-body elastic S-matrix 
of the LL model is [H, [l| 

where p is the momentum difference of the two particles. Consider now the Sh-Gordon model in (1+1) 
dimensions, i.e. the integrable and relativistic invariant field theory defined by the Lagrangian 



^ShG — 2 



2 

2 



^ cosh(g</.) , (18) 

r 



where </> = (^{x, t) is a real scalar field, is a mass scale and c is the speed of light. The parameter /x is 
related to the physical (renormalized) mass m by /^^ = Tram^c^/fi^ sin(7ra), where a — he / {8tt + he g'^) 
pi ]. Let's express the dispersion relation of the energy E and the momentum P of a particle as i? = 
mc^ coshfl, P = mc sinhfl, where 6 is the rapidity and c the light speed. The Sh-Gordon is an integrable 
model and therefore all its scattering processes are purely elastic and can be factorized in terms of the 
two-body S-matrices [l^ 

,^ , sinh — z sin(Q;7r) 
smht^ + I sm(a7r) 

where 9 is the rapidity difference of the two particles. The key observation is that taking simultaneously 
the non-relativistic and weak-coupling limits of the ShG model such that 

g->0, c->oo, p c = 4\/A//i = fixed , (20) 

its 5-matrix (|19p becomes identical to the ^-matrix ([T7)) of the LL model. In this double limit procedure 
the coupling A does not need to be small, i.e. with this mapping we can study the LL model at arbitrarily 
large values of the dimensionless coupling 7. The mapping between the S'-matrices of the two models 
suggests that this mapping should extend also to their Lagrangians and Thermodynamic Bethe Ansatz 
(TBA) equations. This is indeed true. In the non-relativistic limit, the real scalar field can be expressed 
as 



and, once inserted into the Lagrangian (|18p . in the limit c 00 one omits all oscillating terms. More- 
over, when 5 — !■ of eqn (j20[) is considered, the "0^^ terms coming from the potential and kinetic 
parts of cancel each other, while all higher terms of the series expansion of the potential, but the 
quartic one, vanish. Hence, the ShG Lagrangian ()18p reduces to the non-linear Schrodingcr Lagrangian 
P^ . The commutation relation [(l){x,t),Il{x' ,t)] = ihS{x — x') implies for the non-relativistic operators 
mx,t),^Hx',t)]=S{x-x'). 

Similarly the TBA equations of the ShG model (given for instance in [4y|) reduce to the ones of the 
LL model, written down in [l5|. In the LL model at a finite T, the pseudo-energy e{T,ii) satisfies the 
non-linear integral equation 

e{T, M) = ^!Z|!!i^ - ^ o log (1 + e-) , (21) 

where /i is the chemical potential conjugated to the finite density n of the gas, f{p) = ^og Sll{p) is 

the derivative of the phase shift and (fio f = ^ Lpfj) — p')f[p'). The solution of this integral equation 
leads to the free energy and to all other thermodynamical data of the model. 

In the light of the mapping between the S'-matrix, the TBA and the operators of the LL and the Sh- 
Gordon models, we can now proceed to compute the expectation values of the LL model. At equilibrium 
the expectation value of an operator O — 0{x) at temperature T and at finite density is given by 

Tr(e-(^-^^'/('=^'^)0) 
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FIG. 3: gi and 172 using form factors up to 71 = I, 6 and 8 particles, respectively with green dot-dashed, blue 
dashed and red dotted lines. The exact values are given by the solid lines whereas the dot-dashed line below, 
indicated by the arrow, corresponds to the strong coupling expansion ([29}. 

In a relativistic integrable model the above quantity can be neatly expressed as (47j 



(°>^E;^/(ntm)) 



;'^P(0)l^>co„n, (23) 



where f{Oi) = l/(l-|-e^*^^'^) and 9 = 9i, . . . ,6n {0 = dn, ■ ■ ■ , &i) denote the asymptotic states entering the 
traces in ([22]). This formula employs both the pseudo-energy e{9) and the connected diagonal form factor 

of the operator O, defined as (V|e'|'^)^„„„ = FP (^lim^^_^o (OlOl"^, V - ^tt + iV)) where V = Vn, ■ ■ ■ ,m 
and FP in front of the expression means taking its finite part, i.e. omitting all the terms of the form 
rii/rjj and 1/rfl where p is a positive integer. For a recent proof of this formula see (33 |. 

To compute the expectation values of the LL model wc have then to apply eqn ([25)) according to the 
following steps: (a) solve the integral equation (|2ip for e{9)] (b) identify the relevant form- factors of the 
ShG model; (c) take the non- relativistic limit of both the form- factors and eqn ([23| . To proceed, let's 
assume we have solved numerically eqn ([2T|) (a task that can be easily done) and therefore let's attack 
the last two points. The generic m-particle form factor of a local operator O in the ShG model can be 
written as j44 . l45j 

F°(0i, . . . , 0„) ^QZi^u.-., xm) n ^^^^ . (24) 

i<j ■> 

where Xi = e^* and Qjri ^'^s the symmetric polynomials in the x's that fully characterize the operator O. 
The explicit expression of i^min(^) is given in [4J] but the only thing needed here is its functional equation 

„ / „^ „ /„N sinh0 

Fmi„(*^ + e)^^mi„(e) - ■ - ■ I . ■ 

smli 9 + 1 sin(7ra j 

We are interested, in particular, in the symmetric polynomials Qm^ of the exponentials Oq = e'^^'^ since, 
using their Taylor expansion, we can extract the form factors of all normal ordered operators :(/)'':. Their 
close expression is [45{ 



4sin(7rQ;) 

77 



- [g] ( ^^T'"' ) detAUq), (25) 



where M„i{(i) is an (m — 1) x [ra — 1) matrix with elements [Af,„(g)]^^- = (y'^ljli — j + q\- Above, 

\x\ = sin(x7ra)/ sin(7ra) while cri'"'' (a = 0,1,..., m) are the elementary symmetric polynomials in m 
variables. 

We can now compute the local fc-particle correlation function gk of the LL model defined by 

(^tfc^fe) (26) 



where k is an integer (fc = 1,2,3,... ). The ^^'s are functions of the dimcnsionless LL coupling 7 and of 
the reduced temperature r. The relation between gk in the LL model and the corresponding quantity in 
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FIG. 4: g2 at T = 1, 10 using form factors up to ti = i, 6 and 8 particles with green dot-dashed, blue dashed and 
red dotted hues, respectively. The solid hues show the exact result. 



the ShG model in the limit (pO)) is given by 



2m) \ k 



Using eqn f|23p and the connected form factors of the corresponding operator the final expression is 



where 



^2tcon„(te}) = _ , J :^V<co„„({^. = ^}) 



oo,g-i.o ymc / ' mc 

are the double limit (PU)) of the connected form factors. It is worth pointing out that the matrix elements 
of the LL models computed by the Algebraic Bethe Ansatz [s^ can be put in direct correspondence with 
the Form Factors of the relativistic Sh-Gordon model [s^. As shown in the series ([27|) are nicely 
saturated by the first few terms for sufficiently large values of 7 (7 = is a singular point of the model 
[3], therefore one cannot expect a priori any fast convergence nearby). A first check of eqn (|27p is given 
by the case fc = 1: using (|27p (with a chemical potential /i that ensures the finite density n) and summing 
up the series, one easily checks that (ip^ip) = ^nd gi = 1 (a result which comes from translational 
invariance). As shown in Fig. [3l the exact value gi — 1 (solid line) is rapidly approached by just the first 
terms of ^7} : the convergence of the series is always remarkably fast for all 7 > 1.5, where the exact value 
is obtained within a 5% accuracy just using its first four terms. Eq. \27\i also permits to easily recover the 
leading order of the strong coupling (i.e. large 7) expansion of all gk- this always comes from the first 
non-zero integral in the series (|27p and therefore 

^1 /^\l^{k-l) 



9>^ = ^[-) h + (28) 

where Ik = J^-^ dki . . . J^^ dkk Y[i<ji^i ~ ^jY ■ This result coincides with the one obtained in [4lj. 

The quantity g2 can be exactly determined via the Hellmann-Feynman theorem [4^ : its plot at T = 
is shown in Fig. [3] together with our determination from eqn ([27|) and, as before, one observes also in 
this case a fast convergent behaviour of the series. The strong coupling regime of g2 can be computed by 
expanding in powers of 7"^ all the terms in eqn (j27p and for T = we get 

32 = fl - - + (24 - + 0{j~') , (29) 



3 7^ \ 7 57 

in agreement with the result of the Hellmann-Feynman theorem [4l|, \^ . Expression ([29)) is also plotted 
in Fig. [3] in order to show that the determination of 92 (at finite 7) obtained from the first terms of eqn 
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FIG. 5: gz vs 7 at : (a) r = 0, (b) t = 1 and (c) r = 10. In (d) we plot 33 vs r at 7 = 7. In all figures the blue 
dashed and the red dotted lines refer to 71 = 6 and 8 particles, respectively. The solid line in (a) is the exact value 
of (?3 at T = 0. 



(|27p is closer to the exact result, because any of them contains infinitely many powers of 7. At finite 
temperatures the convergence of the series is also pretty good and the results are shown in Fig. |4l 

A new important quantity get by the method above is 33, a quantity known exactly at T = (43l |. 
but only approximately at T > |4l| . This quantity is related to the recombination rate of the quasi- 
condcnsatc and thus to the lifetime of the experiments. Its strong coupling limit at T = is 

16 TT^ / 10 

1576 V 7 



9^ = --[^--]+0{l-^). (30) 



The plot of (73 at r = using form factors up to n = 6 and 8 particles (i.e. one or two terms of the series 
(|27)) ) is in Fig. [^a) and, as in previous examples, it shows a nice convergent pattern to the exact value 
found in [i^. Figs.[5ljb,c) show 53 as a function of 7 at fixed temperature r, while Fig.[5ljd) shows instead 
(73 as a function of r at a fixed value of 7. 

In closing this section, it should be mentioned that the method discussed above works also in the 
case of the super Tonks-Girardeau gas [l^, i.e. a highly excited, strongly correlated state obtained in 
quasi one-dimensional Bose gases by tuning the scattering length to large negative values of 7 by using 
a confinement-induced resonance. 



III. NON-INTEGRABLE ASPECTS 



Two-dimensional massive Integrable Quantum Field Theories (IQFTs) have proven to be one of the 
most successful topics of relativistic field theory: in addition to the examples discussed in the previous 
section, it has found a large variety of applications to many other models (see, for instance, [20|). However, 
despite the elegance and the undeniable success of these methods, the generic situation that occurs in 
many-body quantum physics is that of a non-integrable dynamics: many interesting models fall within 
this class and therefore it would be highly desirable to develop an appropriate formalism to deal with the 
lack of integrability. The breaking of integrability is expected to considerably increase the difficulties of 
the mathematical analysis, since the scattering processes of the excitations are no longer elastic. Non- 
integrable field theories are in fact generally characterized by particle production amplitudes, resonance 
states and, correspondingly, decay events. All these features strongly effect the analytic structure of the 
scattering amplitudes, introducing a rich pattern of branch cut singularities, in addition to the pole 
structure associated to bound and resonance states. 

Among the many possible ways in which a model can be non-integrable, a particularly interesting 
situation occurs when we consider Hamiltonian defined in terms of a conformal Hamiltonian deformed 
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by two relevant operators, eaeh of them giving rise individually to an integrable model 

H = HcFT + ^i / dxipi{x)+X2 / dxtp2{x) . (31) 



There are in fact many interesting physical systems that belong to this class of non-integrable models. 
Let's briefly discuss two of them. 

• The first is the Ising model at A 7^ 1 in an external magnetic field h. We have previously seen 
that when h = 0, the above action corresponds to the integrable theory of the thermal deformation, 
that has only one particle excitation and an elastic 5- matrix equal to 5* ~ —1. On the contrary, 
when T ^ Tc one recovers the integrable theory of the Ising in a magnetic field, whose spectrum 
consists of the 8 massive particles, eq. (jlOp . One may wonder how the spectrum evolves by moving 
the couplings, a question whose answer will be given in the next section. 



The second example is provided by the multi-frequency Sine- Gordon model, with Hamiltonian 
density 

= ^ [{dtV? + {dxvf] + Ai cos (3ip + A2 cos . (32) 

When Ai = 0, this Hamiltonian gives rise to the integrable theory of the Sine-Gordon model with 
frequency a. In addition to the soliton states, such a theory has a number of neutral bound states 

where = ^/(l — a^/Sir). Viceversa, if A2 — 0, we have again a Sine- 



given by[62[ 7V2 



?/3 



When 



Gordon model but of frequency /? and a different number of neutral bound states, Ni 

the ratio of a and /? is a rational number, the potential of the theory has an infinite number of 
periodic and degenerate vacua. On the contrary, when the ratio of the frequencies is an irrational 
number, the potential has only one vacuum that can always be placed at the origin. Also in this 
case, it is interesting to determine how the spectrum evolves by moving the two couplings. 

Notice that, to study the theories associated to an Hamiltonian as (|31|) . it is convenient to regard it as a 
deformation of an integrable Hamiltonina rather than as a multiple deformation of a conformal theory. 
By taking this point of view and grouping differently the terms, the Hamiltonian (|3ip can be written as 

H = H^^, + X, [ dx^,{x) . (33) 



(i = 1, 2, j 7^ i). There are several advantages in doing so. 

1. The first convenience is that the non-integrable theory can be analyzed starting from the basis 
of the Hilbert space provided by the particle excitations associated to the integrable model -fffnt- 
Although the spectra of H and i?'„( would be different, the basis provided by the particles of the 
integrable model will be certainly more appropriate than the conformal one, as far as the infrared 
properties of the non-integrable model are concerned. 

2. The second advantage consists of the exactly solvability of the integrable models, in particular, 
the possibility of computing exactly all the matrix elements (Form Factors) of local and non-local 
operators of such theories. Hence, in a complete analogy with ordinary Quantum Mechanics, one 
can conceive to set up a perturbative approach based on the Form Factors of the integrable models. 
As we will see in the next Section, this perturbative approach will enable us to reach a remarkable 
series of predictions about the mass correction and the decay processes. 

3. When each deformation is individually integrable, there is the obvious freedom of using any of them 
as a starting point. By doing this choice we select a particular basis of the particles and bound states 
thereof. However, since the actual dynamics of the model should be insensitive to such a choice, 
there should be a series of mathematical identities that links one perturbative series to the other. 

In the following we will discuss two different approaches to study the breaking of integrability: the Form 
Factor Perturbation Theory, the Semiclassical Method and the Truncated Conformal Space Approach. 
Each of them provides different information on non-integrable aspects. 
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IV. FORM FACTOR PERTURBATION THEORY 



For massive non-integrable field tlieories, a convenient perturbative scheme was originally proposed in 
@ and called Form Factor Perturbation Theory (FFPT), since it is based on the knowledge of the 
exact Form Factors of the original integrable theory. Its generalization to massless case has been discussed 
in [s^l. In these papers it was shown that, even using just the first order correction of the FFPT, a great 
deal of information can be obtained, such as the evolution of their particle content, the variation of their 
masses and the change of the ground state energy. Whenever possible, universal ratios can be computed 
and successfully compared with their value obtained by other means, as the universal ratios relative to 
the decay of the particles with higher masses in the Ising model in a magnetic field, once the temperature 
is displayed away from the critical value 1321 . For other and important aspects of the Ising model along 
non- integrable lines see the references [3, ESmOl l . 



A. First Order Perturbation Theory 



Referring the reader to the original papers 0, Q for the theory of Form Factor Perturbation Theory, 
here we report the conclusions which can be reached by studying the first order correction to theory as 
the one of cq. ([31]) . where the integrable model is defined by the action ((3T|) with A2 = 0, whereas its 
perturbation is given by the relevant scalar operator (p2{x). Let xi and X2 be the scaling dimensions of the 
two operators ipi and (p2. The theory depends in this case on the two dimcnsionful couplings constants 
Ai and A2 Since Ai ^ M^^^^ and A2 ^ M^^^^ (where M is a mass scale), we can decide to use Ai 
as dimcnsionful parameter of the theory and the dimensionless combination 

X = X2xf^ (34) 

as a label of the different Renormalization Group trajectories which originate from the fixed point at 
Ai = A2 = 0. For example, if N{x) denotes the number of stable particles in the spectrum of the theory, 
their masses can be expressed as 

ma{Xi,x),= Ca{x)xf^ , a = l,2,...,iV(x) , (35) 

where Ca{x) is an amplitude which characterises the whole trajectory. Similarly, the vacuum energy 
density can be written as 

fvac(Ai,x)=f(x)Ar^ . (36) 

Dimensionless quantities, as for instance mass ratios, only depend on x a-nd therefore they do not vary 
along the trajectories of the Renormalization Group. 

Once the new interaction A2 j cP'x(p2{x) is switched on in the action, the integrability of the unper- 
turbed theory is generally lost and the S-matrix amplitudes become complicated quantities. Inelastic 
processes of particle production are no longer forbidden and, as a consequence, the analytic structure of 
the scattering amplitudes present additional cuts due to the higher thresholds. In particular, their ex- 
pression is no longer factorized into the sequence of two-body scattering amplitudes and, even in elastic 
channels, the only surviving restriction on the final momenta comes from energy-momentum conservation. 

The knowledge of the matrix elements of the perturbing field 'f2{x) ensures the possibility to compute 
perturbatively both the amplitudes of the inelastic processes and the corrections to the elastic ones. For 
instance, the first order corrections to the masses of the particles and to the vacuum energy density are 
given by 

5Ml 2A2 i^^; [i^, 0) (5™„^, , (37) 



Jfvac ^ A2 [o(0|</52|0)o] . (38) 

The best use of these formulas is to get rid of the explicit dependence on the normalisation of the 
perturbing operator by defining universal quantities, as ratios of the mass shifts. Hence, under the validity 
of the linear approximation, all the universal quantities of non-integrablc field theories can be entirely 
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expressed in terms of Form Factors of the integrable ones. A comparison of the theoretical predictions with 
their numerical determinations will be presented in the sequel, when we discuss the class of universality 
of the Ising model. Let's discuss now the important link which exists between the confinement of the 
excitations and the non-local nature of the perturbing operator. 

Non-locality and Confinement Let's consider in more details the mass correction given in eq. ([37)) 
where the the form factor of the operator 'f2{x) is defined by the matrix element 

F:i{e)^{Q\^2{Q)H0M02)) ■ (39) 
Let us recall that the 2-particle Form Factor of an integrable theory satisfies the equations 

Fg{e) = sfMFg{-e) , (40) 



In the second equation the explicit phase factor e^'^'^^"'"-'^ is inserted to take into account a possible semi- 
locality of the operator which interpolates the particles and the operator 0{x). If 7 7^ 0, the 2-particle 
Form Factor presents a pole at 9 = ±i7r, with the residue given by 

~ ^Rese=±^^F^M - (1 - e^^''^^o-^){0\O\0) . (42) 

Hence, according whether the perturbing field is local or non-local with respect to the asymptotic parti- 
cles, there are two different scenarios. If the field that breaks intcgrability is a local operator, the mass 
correction of the particles is finite. Viceversa, if the perturbing field is non-local, the mass correction of 
the particles is divergent. The last case implies the confinement of the particles, that occurs as soon as 
the non-integrable perturbation is turned on. 

There are several ways to show the confinement phenomena. One consists of computing the propagator 
{A(jj)A(—p)) of the particle A in the perturbed theory. At the tree level approximation, this consists of a 
geometrical scries that can be explicitly summed and for the propagator of the perturbed theory we have 

If Srn^ — 00, the propagator obviously vanishes, i.e. the particle disappears from the spectrum. A more 
intuitive explanation of the confinement phenomenon comes from the analysis of the Ising model, discussed 
below. 

Let's close this section by mentioning that, applied to the double Sine-Gordon model @, the FFPT has 
helped to clarify the rich dynamics of this non-integrable model. In particular, in relating the confinement 
of the kinks in the deformed theory to the non-locality properties of exponential operators with respect 
to the the kinks and in predicting the existence of a Ising-like phase transition for particular ratios of the 
two frequencies - results which were later confirmed by a numerical study [sij . The FFPT has been also 
used to study the spectrum of the 0(3) non-linear sigma model with a topological 9 term, by varying 9 
[sgI \3T\ , a model closely related to anti- ferromagnetic spin chains f38l - l40| . 



B. The scaling Region of the Ising Model 



It is interesting to study the evolution of the mass spectrum of the Ising model by moving its couplings 
along the path C shown in Fig. [5] in the plane (r, /i) [5[: this curve starts from the low temperature 
phase of the model and ends at its high-temperature phase, here represented by the points (a) and (f) 
respectively. The action of the model is given by 

A = AcFT + r j (fxe{x)+h j (fxa{x) , (44) 

and gives rise to a massive integrable model when one of the two coupling constants is switched off. The 
action defines a family of field theories identified by x = t|/i|~^/^^ G (—00,-^00), a dimcnsionless 
RG invariant quantity. The spectrum of the theory changes in a significant way by moving x- 
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h 




FIG. 6: Interpolation curve in the plane (r, h) between the low and high-temperature phases of the Ising model. F 
is a Renormalization Group trajectory, identified by the dimensionless parameter X- 



In the low-temperature phase (corresponding to x = ~c« and to the point (a) of the curve C), the 
model has two degenerate vacua and therefore its excitations consists of the topological kink and anti- 
kink that interpolated between the two ground states. Along the magnetic axes (x = 0, corresponding 
to the point (d) of the curve C), the spectrum of the model consists instead of 8 particles, with different 
masses. Finally, in the high-temperature phase (i.e. x = the system has a unique vacuum and only a 

massive excitation above it. Let's see how this scenario can be recovered by the Form Factor Perturbation 
Theory. 

Let's start our analysis from the point (a), where the massive excitations are the kink/anti-kinks that 
interpolate between the two degenerate vacua. By switching on the magnetic field, the model moves to 
the point (b) of the curve C. The form factor of the perturbing field a on the two-particle kink/ ant i- kink 
state is given by 

F^0i-O2) = (O|a(O)|A(0i)A(02)) = tanh^^y^ . 

Therefore eqn. ()37|) leads to an infinite correction to the mass of the kinks, i.e. the kinks get confined 
as soon as the magnetic field is switched on. Looking at the effective potential of the theory, it is not 
difficult to see that this is the correct conclusion: no matter how small the magnetic field may be, it lifts 
the degeneracy of the two vacua, as shown in Fig. [3 and consequently there is no longer the possibility 
of having topological configurations. 

Consider now the effect of the magnetic field on a state made of a kink and an antikink separated by 
a distance R. When the magnetic field is absent, the energy of this state is essentially equal to 2M, i.e. 
the sum of the masses of the kink and the antikink. The energy of this state depends very weakly on 
the distance R because this field configuration takes values on the zeros of the effective potential and, no 
matter how large the distance R could be, there is no change in the energy of this state. This situation 
changes by switching on the magnetic field since, in this case, at every point of the space there is an energy 
gap equal to 2h and the energy U of this state becomes a linear function of R, U{R) = 2AI + 2hR. This 
attractive interaction between the kink and the antikink gives rise to a discrete spectrum of bound states. 
Regarding the kinks as very massive and quasi-static, the energy of the bound states can be obtained 




(a) (b) 

FIG. 7: Ejfecttive potential m the low-temperature phase (a) and in the presence of an infinitesimal magnetic field 
(h). In the last case, the two minima are no longer degenerate and the kink/antikink disappear from the spectrum 
of the asymptotic states. 
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by solving the quantum mechanical problem of the bound states for a linear potential, well-known in 
Quantum Mechanics. The result is simply 

Eu = mu - (2 + h^'^-il''')M , (45) 
where 7fc are the positive roots of the equation 

Jilk) = Q7fc) + J- 1 Qt^) = , 

{Jy{x) is the Bessel function of order v). The structure of the bound states is shown in Fig. [51 Obviously 
not all these states are stable: the stable ones are identified by the condition m„ < 2mi, while all particles 
with a mass higher than the threshold 2mi decay into particles of lower masses. When x increases, i.e. 
when the system moves clockwise along the curve C of Fig. [51 the number of bound states monotonically 
decreases. At the point (d), there are the 8 stable particles of the Zamolodchikov's solution of the Ising 
model in a magnetic field. 

It is worth stressing that in this case the five particles with mass higher than the threshold are stable 
just for the integrability of the model. Moving away from the magnetic axes by means of the operator 
e{x), the first three particles change the value of their masses, while the remaining five particles decay 
into the low-energy channels. To estimate both effects, we need the Form Factors of the energy operator 
but in the integrable theory of the Ising model in a magnetic field, determined in Q. Here we simply 
report their expressions 

(0|e(0)|0) = mi , 

Ffi(f7r) = {{)\e{{))\Ai{e + iT:)Ai{d)) = -17.8933. .mi , 
F^^tt) = {Q\e{Q)\A2{e + i-K)A2{e)) = -24.9467..mi , 
F|3(i7r) = {Q\e{Qi)\A3{e + m)A3{e)) = -53.6799. .mi , 



The first equation may be regarded as the normalization condition of the energy operator e{x). The 
corrections of the universal ratios are given by 





(0|e|0) 


m° 


Sm2 




mi 


Smi 




m§ 


5mz 




mi 


Smi 




ml 



-0.0558..m; 



— i = 0.8616.. 



= 1.5082.. 



(46) 



In turn, these quantities can be independently determined by a numerical solution of the model [11 and 




(a) 



(b) 



FIG. 8: (a): kink-antikink state separated by a distance R; (h) kink- antikink potential in the presence of a magnetic 
field h and its hound states. The stable hound states are identified by the condition En < Et- 
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abc \) abc 




FIG. 9: Perturbative diagram at the second order in r relative to the imaginary part of the mass of the particle 
c. The intermediate particles a and b satisfy the on-shell conditions p\ = m\ e pi — ml. When c > 5, there are 
additional diagrams with more intermediate particles. 



the values determined in this way are 

5£,, 



Smi 
Sm2 
Smi 
5mz 



- -0.05 , 

0.87 , (47) 
1.50 . 



8mi 

As it can be seen from the expressions above, there is a satisfactory agreement between the theoretical 
and the numerical estimates. 

Breaking the integrability of the Ising model in a magnetic field has a more dramatic effect on the five 
particles with a mass above threshold. Their stability is only due to integrability and, in its absence, they 
decay. These decay processes were studied in [32| . In the perturbative approach, the decay processes are 
associated to the presence of a negative imaginary part in the mass that is a second order perturbative 
effect in r, as shown in Fig. HI 



2 _ ™ -n _ ^, _2 ol-<5a6 I /cafe P 

sinh^i^""' 



Imm^ = - ^ iricTc^ab — -t'^ ^ 



where Ti^^ab is the decay amplitude of the particle A^. into the two particles Aa A;,, whereas 

Uab = F!^,{^^,eir')Mr'^) 

T — 

The rapidities ^i^"^^ and 0^^°'''^ are fixed by the conservation of the energy and the momentum in the 
decay process Ac — > AaA^ in the rest frame of the particle Ac. In the above equations all masses are 
the unperturbed values at t = 0. When c > 5 the sum must be completed including the contribution of 
the decay channels with more than two particles in the final state. Once the decay amplitudes Vc^ab are 
known, one can determine the lifetime tc of the unstable particle Ac given by 

ic-^, r, = ^r,^,b. (48) 

a<b 

For the Ising model, the relevant matrix elements are 

|/4n| = (36.73044..) |(e)|^^o 
l/siiH (19.16275..) |(e)|,^o 
I/512H (11.2183..) . 

where the normalization of the operator e is fixed by its vacuum expectation value 

{e)r=o = (2.00314..) \h\'^'' . 

The imaginary part of the mass of the first two particles, which are over threshold, is given by 

Imm^ ~ (-840.172..) fli£ll^\ = (-173.747..) 



nil 



Immi ~ (-240.918..) ( lifZl^ ) = (-49.8217..) 
* nil 



2 
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The ratio of their hfetime is universal 

li,, t4 ^ j.^^ m4lmm| ^ ^ 23326.. (49) 

While the particle A4 can only decay into AiAi, the particle A5 can also decay into the channel ^1^2- 
The ratios of the amplitudes of these decays 

TOc|r=0 Tc-s-afc 

Oc^afc — 



are given by 



llm m?| 



limfeg^ii = 0.47364.., lim 65^12 0.52635. 

r— 5-0 r-i-O 



Notice that eq. (HH]) predicts that the lifetime of the particle A^ is almost four time longer than the 
lifetime of the particle A^. This paradoxical result, in contradiction with the intuitive idea that a heavy 
particle should decay faster than a light one, finds its explanation once again in the peculiar behavior of 
the phase space in two dimensions. For the decay process Ac — AaAb the phase space in d-dimcnsion is 
given by 

(50) 

Pa PI "Ic 

where p ~ \\p)^\ ~ \\p)b\ value in the rest frame. For fixed decay products, p grows with 771^: in 

d = 2, this term joins the factor nic in the denominator and leads to a suppression of the phase space. In 
the Ising model, cqns. (|49 p - (|49p show that this suppression is further enhanced by the dynamics (i.e. by 
the vaues of the matrix elements) in a way that is not compensated by the additional decay channels. 

If we keep moving along the curve C, one firstly meets a value xi a-t which the mass of the particle 
A3 becomes larger than 2mi and, later on, a second value X2 at which also the mass of the particle A2 
becomes larger than 2toi. When x > X2 the spectrum of the stable particles of the theory consists of 
only one excitation. In the limit x ~^ +00, this is nothing but the particle of the integrable theory of the 
high-temperature phase of the model. 



V. SEMICLASSICAL METHODS 



Non-integrable quantum field theories with topological excitations (kink-like) can be also studied by 
semiclassical methods. Such theories are described by a scalar real field (p(x), with a Lagrangian density 

= lid^^r - Uiip) , (51) 

where the potential U{ip) possesses several degenerate minima at (/si"' (a = 1, 2, . . . , n), as the one shown 
in Figure 1101 These minima correspond to the different vacua | a ) of the associate quantum field theory. 
Excitations of the system. The basic excitations of this kind of models arc kinks and anti-kinks, 
i.e. topological configurations which interpolate between two neighbouring vacua. Semiclassically they 
correspond to the static solutions of the equation of motion, i.e. 

dl^ix) = U'Mx)] , (52) 




FIG. 10: Potential U{ip) of a quantum field theory with kink excitations. 
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FIG. 11: Kink and antikink configurations. 



with boundary conditions (^(— oo) = ip^^ and ip{+oo) ~ where b 
solutions of this equation, their classical energy density is given by 



1 / dlfab 

2\dx 



U{ipab{x)) 



a ± 1. Denoting by ipab{x) the 



(53) 



and its integral provides the classical expression of the kink masses 

oo 

Mat = / f^abix) ■ 



(54) 



It is easy to show that the classical masses of the kinks (pab{x) are simply proportional to the heights of 
the potential between the two minima (^i"^ and ip''^^ . Once put in motion by a a Lorentz transformation, 
i.e. (pab{x) — )■ <Pab [{x ± vt) / \/ 1 ~ v'^ \ , thcse configurations describe in the quantum theory the kink states 
I Kab{0) ), where a and h are the indices of the initial and final vacuum, respectively. The quantity 9 is 
the rapidity variable which parameterises the relativistic dispersion relation of these excitations, i.e. 



E = Mnh coshf 



P ^Mnh sinh6' 



(55) 



Conventionally | Ka,a+i{9) ) denotes the kink between the pair of vacua {| a), | a + 1 )} while | Ka+i^a ) 
is the corresponding anti-kink. For the kink configurations one can adopt the simplified graphical form 
shown in Figure [TT] The multi-particle states are given by a string of these excitations, with the adjacency 
condition of the consecutive indices for the continuity of the field configuration 



K 



(a 



i+1 



±1) 



(56) 



In addition to the kinks, in the quantum theory there may exist other excitations, i.e. ordinary scalar 
particles (breathers). These are the neutral excitations | Bc{9) )a (c = 1,2,.. .) around each of the vacua 
I a). For a theory based on a Lagrangian of a single real field, these states are all non-degenerate: in fact, 
there are no extra quantities which commute with the Hamiltonian and that can give rise to a multiplicity 
of them. The neutral particles must be identified as the bound states of the kink-antikink configurations 
that start and end at the same vacuum | a), i.e. | Kab{9i) Ki,a{92)), with the "tooth" shapes shown 
in Figure 1121 If such two- kink states have a pole at an imaginary value i u'^^ within the physical strip 
< ImO < TT of their rapidity difference 9 = 9i — 02, then their bound states are defined through the 
factorization formula which holds in the vicinity of this singularity 



Kabi9l)Kbai92)) 



9ab 



Br. 



(57) 



io> 



io> 



I0> 



FIG. 12: Kink-antikink configurations which may give rise to a bound state nearby the vacuum \ 0)a. 
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In this expression g^f^ is the on-shell 3-particle couphng between the kinks and the neutral particle. 
Moreover, the mass of the bound states is simply obtained by substituing the resonance value i ujjj, 
within the expression of the Mandelstam variable s of the two-kink channel 

s = 4M2j, cosh" ^ ^ mc = 2MabC0S ^ . (58) 

Concerning the vacuum states, in the infinite volume their classical degeneracy is removed by selecting 
one of them, say | fc ), out of the n available. This happens through a "spontaneously symmetry breaking 
mechanism" , even though there may be no internal symmetry to break at all, as for the potential of Figure 
1101 In the absence of a symmetry which connects the various vacua, the excitations above them may be 
very different and it is interesting to address the question of their and the value of their masses. This 
problem can be solved through the application of a remarkably simple formula due to Goldstone-Jackiw 
[531 : this applies in the semiclassical approximation, i.e. when the coupling constant goes to zero and the 
mass of the kinks becomes correspondingly very large with respect to any other mass scale. In its refined 
version, given in [53 |. this formula reads as follows 

oo 

f^,ie) = {Kabiei) I ^{0) I Kab{e2)) ^ j dxe''--'- ^ab{x) , (59) 

— oo 

where 9 = 9i—62. Substituting in this formula 9 — > in— 6, the corresponding expression may be interpreted 
as the following Form Factor 

Fabi(^) = fi^^-d) = (« I m I Kba{02)) , (60) 

where appears the neutral kink states around the vacuum | a) of interest. Eq. (|59p deserves several 
comments. 

1 . The appealing aspect of the formula ([5^ stays in the relation between the Fourier transform of the 
classical configuration of the kink, - i.e. the solution (pab{x) of the differential equation (|52|) - to 
the quantum matrix element of the field (p(0) between the vacuum | a ) and the 2-particlc kink state 

I Kab{9l)Kbai92)). 

Given the solution of eq. (|52p and its Fourier transform, the poles of Fab{9) within the physical 
strip of 6 identify the neutral bound states which couple to (p. The mass of the neutral particles 
can be extracted by using eq. (|58p , while the on-shcU 3-particle coupling can be obtained from 
the residue at these poles (Figura [T3)) 

hm^ {6 - *<,) FabiO) = igU {a \ ^(0) \ B,) . (61) 



2. For a generic theory, the classical kink configuration ipabix) is not related in a simple way to the 
anti-kink configuration tpba{x). This is the reason why neighbouring vacua may have a different 
spectrum of neutral excitations, as shown in the examples discussed below. 
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3. This procedure for extracting the bound states masses permits in many cases to avoid the semi- 
classical quantization of the breather solutions [52j . making their derivation much simpler. The 
reason is that, the classical breather configurations depend also on time and have, in general, a 
more complicated mathematical structure than the kink ones. Moreover, it can be shown that in 
non-integrable theories these configurations do not exist as exact solutions of the partial differential 
equations of the field theory. On the contrary, to apply eq. (|59p , one simply needs the solution of 
the ordinary differential equation (|52p . It is worth notice that, to locate the poles of f^ij{0), one 
has simply to look at the exponential behavior of the classical solutions at a; — )• ±03, as discussed 
below. 

Below we will discuss a class of theories with only two vacua, which can be either symmetric or 
asymmetric ones. A complete analysis of other potentials can be found in the original paper [ssj . 



A. Symmetric wells 

A prototype example of a potential with two symmetric wells is the ip"^ theory in its broken phase. The 
potential is given in this case by 

Let us denote with | ±1) the vacua corresponding to the classical minima (p^^^^ = ±-^- By expanding 



around them, (p — (p^^^ + rj, we have 



C/(</j^°' + Tj) = T]"^ ± mVXif + ^ri'^ . (63) 

Hence, perturbation theory predicts the existence of a neutral particle for each of the two vacua, with a 
bare mass given by mf, = \/2m^ irrespectively of the value of the coupling A. Let's see, instead, what is 
the result of the semiclassical analysis. The kink solutions are given in this case by 



V~a.a\.x) = " ^^^^ 



V2 



a = ±1 (64) 



and their classical mass is 



Mo = j e{x)dx = • (65) 

—00 

The value of the potential at the origin, which gives the height of the barrier between the two vacua, can 
be expressed as 

UiO)^^Mo, (66) 

and is proportional to the classical mass of the kink. If we take into account the contribution of the small 
oscillations around the classical static configurations, the kink mass gets corrected as (52| 

M^'-^"4^m('')+OiX) . (67) 



Defining 



and the adimensional quantities 



27r 4^3 



7rV2 2V6 



> 
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the mass of the kink can be expressed as 



M 



y/2m 



(69) 



Since the kink and the anti-kink solutions are equal functions (up to a sign), their Fourier transforms 
have the same poles and therefore the spectrum of the neutral particles will be the same on both vacua, 
in agreement with the Z2 symmetry of the model. We have 



dx e 



iMe X 



1 



sinh 



\/2m 



By making now the analitical continuation 9 ^ i 

F-a,a{S) = (a I ViO) I K_a,a{^l)I<a,-a{e2)) « 



9 and using the above definitions (|68p . we arrive to 
1 



sinh 



The poles of the above expression are located at 

9n = in{l - (n) 

and, if 



0,±1,±2,.. 



(70) 



(71) 



^ > 1 , (72) 
none of them is in the physical strip < Im 9 < -k. Consequently, in the range of the coupling constant 

1.02338... (73) 



TO 



A 27r 1 

~ - YTT 



TTC 



the theory does not have any neutral bound states, neither on the vacuum to the right nor on the one to 



the left. Viceversa, if ^ < 1, there are n ■ 



neutral bound states, where [x] denote the integer part of 



the number x. Their semiclassical masses are given by 









sin 




~ n Tfih 







3 2 
^-32;^" 



(74) 



Note that the leading term is given by multiples of the mass of the elementary boson | B\). Therefore 
the n-th breather may be considered as a loosely bound state of n of it, with the binding energy provided 
by the remaining terms of the above expansion. But, for the non-integrability of the theory, all particles 
with mass m„ > 1m\ will eventually decay. It is easy to see that, if there are at most two particles in 
the spectrum, it is always valid the inequality TO2 < 2toi. However, if ^ < i, for the higher particles one 
always has 



TOfc > 2mi 



for /j = 3, 4, 



(75) 



According to the semiclassical analysis, the spectrum of neutral particles of ip'^ theory is then as follows: 
(i) if ^ > 1, there are no neutral particles; (ii) if ^ < ^ < 1, there is one particle; (iii) if i < f < i there 

arc two particles; (iv) if ^ < ^ there are |- particles, although only the first two are stable, because the 
others are resonances. 

Let us now briefly mention some general features of the semiclassical methods, starting from an equiv- 
alent way to derive the Fourier transform of the kink solution. To simplify the notation, let's get rid of 
all possible constants and consider the Fourier transform of the derivative of the kink solution, expressed 
as 



G(fc) 



dx e 



ikx 



1 



cosh^ ; 



(76) 
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FIG. 14: Neutral bound states of (p* 
higher lines are the resonances. 



theory for g < 1. The lowest two lines are the stable particles whereas the 



We split the integral in two terms 



G(fc) 



dx e 



ikx 



cosh 



dx e 



ikx 



cosh^ X 



(77) 



and we use the following series expansion of the integrand, valid on the entire real axis (except the origin) 

1 



cosh^ X 



4 ^(-l)"+ine-2" 



\x\ 



Substituting this expression into (1771) and computing each integral, we have 



Gik) = 4z£(-l) 



1 



1 



ik + 2n —ik + 2n 



(78) 



(79) 



Obviously it coincides with the exact result, G{k) ~ nk/ sinh-|fc, but this derivation permits to easily 
interpret the physical origin of each pole. In fact, changing k to the original variable in the crossed 
channel, k — )• {in — 9)/£,, we see that the poles which determine the bound states at the vacuum | a) are 
only those relative to the exponential behaviour of the kink solution at x — )- — oo. This is precisely the 
point where the classical kink solution takes values on the vacuum | a). In the case of (p"^, the kink and 
the antikink are the same function (up to a minus sign) and therefore they have the same exponential 
approach at x = —oo at both vacua | ±1). Mathematically speaking, this is the reason for the coincidence 
of the bound state spectrum on each of them: this does not necessarily happens in other cases, as we will 
see in the next section, for instance. 

The second comment concerns the behavior of the kink solution near the minima of the potential. In 
the case of (^^, expressing the kink solution as 



(p{x) 



tanh 



V2 



m — 1 



VA + 1 



and expanding around x = — cx). we have 



'Pit) 



Vx 



[1 - 2t + 2<2 - + • • • 2 (-l)"t" • • • ] 



(80) 



(81) 



where t = exp[V2a;]. Hence, all the sub-leading terms are exponential factors, with exponents which are 
multiple of the first one. Is this a general feature of the kink solutions of any theory? It can be proved 
that the answer is indeed positive [55[. 

The fact that the approach to the minimum of the kink solutions is always through multiples of the 
same exponential (when the curvature uj at the minimum is different from zero) implies that the Fourier 
transform of the kink solution has poles regularly spaced by £,a = -^f^r^ the variable 9. If the first of 
them is within the physical strip, the semiclassical mass spectrum derived from the formula (j59|) near the 
vacuum | a ) has therefore the universal form 



m„ = 2Mab sin n 



(82) 
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FIG. 15: Example of (p^ potential with two asymmetric wells and a bound state only on one of them. 



As we have previously discussed, this means that, according to the value of ^a, we can have only the 
following situations at the vacuum | a): (a) no bound state if > 1; (b) one particle if ^ < £,a < 1; 

(c) two particles if ^ < < i; (d) j- particles if < 5, although only the first two are stable, 
the others being resonances. So, semiclassically, each vacuum of the theory cannot have more than two 
stable particles above it. Viceversa, if a; = 0, there are no poles in the Fourier transform of the kink and 
therefore there are no neutral particles near the vacuum \ a). 



B. Asymmetric wells 



A polynomial potential with two asymmetric wells must necessarily employ higher powers than ip'^, 
and the simplest example is given by a polynomial of maximum power ip^. Beside its simplicity, the ip^ 
theory is relevant for the class of universality of the Tricritical Ising Model, and the information available 
on this model turn out to be a nice confirmation of the semiclassical scenario discussed below. 

A class of potentials with two asymmetric wells is given by 

with a,b,c all positive numbers. To simplify the notation, it is convenient to use the dimensionless 
quantities obtained by rescaling the coordinate as — > mx'^ and the field as f^x) — > \/X/mip(x), so that 
the lagrangian of the model becomes 

(84) 

The minima of this potential are localised at ip^^'' = —a and ip'""' ~ b and the corresponding ground states 
will be denoted by | ) and | 1 ) , with the curvature of the potential at these points given by 



m 



1(9^)2 -i(^ + a)2(^- 6)2(^2 +c) 



C/"(-a) ^ a.2 = (a + 6)2(a2 + c) ; 

t/"(6) EE = (a + &)2(62+c) . ^ ' 

For a ^ 6, we have two asymmetric wells, as shown in Figure [151 Let's assume that the curvature at 
the vacuum | 0) is higher than the one at the vacuum | 1), i.e. a > b. The problem we would like to 
examine is whether the spectrum of the neutral particles \ B)s {s = 0,1) may be different at the two 
vacua, in particular, whether it would be possible that one of them (say | 0)) has no neutral excitations, 
whereas the other has just one neutral particle. The ordinary perturbation theory shows that both vacua 
has neutral excitations, although with different value of their mass: 

m(°) ia + b)y/2 (a^ + c) , m^^) = {a + b)y/2 {b^ + c) . (86) 

Let's see, instead, what is the semiclassical scenario. The kink equation is given in this case by 



-f = ±(^ + a)i^ - b) V^T~c . (87) 
ax 

Even in the absence of an exact solution, we can extract the main features of this equation. The kink 
solution interpolates between the values —a (at x ~ —00) and b (at x = +00), while anti-kink solution 
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docs viccversa, but with an important difference: its behaviour at a; = — oo is different from the one of 
the kink. As a matter of fact, the behaviour aX x = — oo of the kink is always equal to the behaviour at 
X = +00 of the anti-kink (and viceversa), but the two vacua are approached, in this theory, differently. 
This is explicitly shown in Figure [1^] and proved in the following. 

Let us consider the limit x — > —00 of the kink solution. For these large values of x, we can approximate 
eq. (|87p by substituting, in the second and in the third term of the right-hand side, (p ~ —a, with the 
result 



i'^] ~ {ip + a){a + b)\/a^ +c , 



This gives rise to the following exponential approach to the vacuum | 0) 

ip^) i(x) ~ —a + Acxp(a;ox) 



X -> 



(89) 



where A > is a arbitrary costant (its actual value can be fixed by properly solving the non-linear 
differential equation). To extract the behavior at a; — >■ —00 of the anti-kink, we substitute this time (ys ~ 6 
into the first and third term of the right hand side of (|87l) . so that 

x^ -00 (90) 



This ends up in the following exponential approach to the vacuum | 1) 

(^i_o(x) ~ 6 — B exp(a;ia;) , 2; — > — 



(91) 



where i? > is another constant. Since luq 7^ wi, the asymptotic behaviour of the two solutions gives rise 
to the following poles in their Fourier transform 



-^(■^0,1) 



A 



ik 



(92) 



-B 



(jji + ik 



In order to locate the pole in 9, we shall reintroduce the correct units. Assuming to have solved the 
differential equation (j87p . the integral of its energy density gives the common mass of the kink and the 
anti-kink. In terms of the constants in front of the Lagrangian (|84p . its value is given by 



AI 



A2 



(93) 



where a is a number (typically of order 1), coming from the integral of the adimensional energy density 
((54)) . Hence, the first polefi^] of the Fourier transform of the kink and the antikink solution are localised 
at 

A2 



0(0) 



ITT [1 — LOq 



ZTT ( 1 — CJl 



m \ 



ttMJ 



ITT [ 1 — UJQ ■ 



in [ 1 — u)i 



A2 



(94) 




FIG. 16: Typical shape of (ff )gj, obtained by a numerical solution of eq. \87^ . 
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If we now choose the couphng constant in the range 

- < 4 < - ' (95) 

Wo TO LOi 

the first pole will be out of the physical sheet whereas the second will still remain inside it! Hence, the 
theory will have only one neutral bound state, localised at the vacuum | 1 ) . This result may be expressed 
by saying that the appearance of a bound state depends on the order in which the topological excitations 
are arranged: an antikink-kink configuration gives rise to a bound state whereas a kink-antikink docs not. 

Finally, notice that the value of the adimensional coupling constant can be chosen so that the mass of 
the bound state around the vacuum | 1 ) becomes equal to mass of the kink. This happens when 

Strange as it may appear, the semiclassical scenario is well confirmed by an explicit example. This 
is provided by the exact scattering theory of the Tricritical Ising Model perturbed by its sub-leading 
magnetization. Firstly discovered through a numerical analysis of the spectrum of this model [56j . its 
exact scattering theory has been discussed later in (57| . 



VI. CONCLUSIONS 



In this paper we have shown that quantum integrablc models have enough set of constraints to allow the 
determination of their mass spectrum and correlation functions. Quantum Ising model and Lieb-Liniger 
model have been our paradigmatic examples to illustrate several aspects of quantum integrability. We 
have discussed how to control the breaking of integrability by means of two approaches: the Form Factor 
Perturbation Theory and the semiclassical method. When the breaking of integrability is realized by 
a non-local operator, one gets the confinement phenomena of the topological excitations. Form Factor 
Perturbation Theory permits also to control the decay processes and computation of the life-time of 
unstable particles. In addition, we have used simple arguments of the semi-classical analysis to investigate 
the spectrum of neutral particles in quantum field theories with kink excitations. We have analyzed, in 
particular, two cases: the first relative to a potential with symmetric wells, the second concerning with a 
potential with asymmetric wells. As a general result of this analysis, we have the existence of a critical 
value of the coupling constant, beyond which there are no neutral bound states. Another result is about 
the maximum number n < 2 of neutral particles above a generica vacuum of a non-integrable theory. 
An additional aspect is the role played by the asymmetric vacua, which may have a different number of 
neutral excitations above them. Let's finally mention that integrable and non-integrable models can be 
numerically studied by means of the so-called Truncated Conformal Space Approach [581 , underlying for 
instance the different statistics which rules the energy level difference, as shown in [59| 
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